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ABSTRACT 



We present the results of a high resolution study with the Submillimeter Array 
towards the massive star forming complex G9. 62+0. 19. Three sub-mm cores are de- 
tected in this region. The masses are 13, 30 and 165 M© for the northern, middle and 
southern dust cores, respectively. Infall motions are found with HCN (4-3) and CS 
(7-6) lines at the middle core (G9.62+0.19 E). The infall rate is 4.3 x 10"^ MQ-yr'^. 
In the southern core, a bipolar-outflow with a total mass about 26 M© and a mass-loss 
rate of 3.6 x 10~^ Mg-yr"^ is revealed in SO (87 - 7?) line wing emission. CS (7-6) 
and HCN (4-3) lines trace higher velocity gas than SO (87 - 77). G9.62+0.19 F is 
confirmed to be the driving source of the outflow. We also analyze the abundances of 
CS, SO and HCN along the redshifted outflow lobes. The mass-velocity diagrams of 
the outflow lobes can be well fitted by a single power law. The evolutionary sequence 
of the cm/mm cores in this region are also analyzed. The results support that UC Hii 
regions have a higher blue excess than their precursors. 

Subject headings: Massive core:pre-main sequence-ISM: molecular-ISM: kinematics and 
dynamics-ISM: jets and outflows-stars: formation 
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1. Introduction 



High-mass stars play a major role in th e evolution of the Galaxy. 



sources of heavy elements and UV radiation dZinnecker & Yorke 



They are the principal 



20071). However, the formation 



and evolution of high-mass stars are still unclear. A possible evolution s equence of high-mass 
stars from infrared dark clouds to classic Hii regions has been suggested dVan der Tak & Menten 
20051). But one of the major topics whether high-mass s tars form through accretion-disk- 



outflow, like low-mass ones (.Shu. Adams & Lizano 



(|Wolfire. & Cassinelli 



1987 



Bonnell et al 



1987|) . or form via collision-coalescence 



19981) is still far from solved. 



Yet more and more observations at various resolutions seem to support the accretion- 
disk-outflow models rather than collision-c oalescence models. Disks are detected in several 
high- mass star forming regions (|Patel et al.l 



2005 



Jiang et al. 



2005 



Sridharan. Williams. & Fuller 



20051). Outflows are found with a high de t ection rate as in low-mass cores in single-dish 



surveys dWu et al. 



2004 



Zhang et al. 



2005 



Oin et al 



2008a|) . High resolution studies have 



also confirmed that molecular outflows are common in high-mass star form ing regions 



(ISu. Zhang. & Lim 



2004 



Oiu et al. 



20071 : 



Oin et al. 



2008b 



c; 



Oiu et al. 



2009h. Search ing 



for inflow motions also has made 



Fuller. Williams. & Sridharan 



2007 



2OO9I : 



2005 



arge progress in recent vears (.Wu & Evans 



Wvrowski et al. 



Furuya. Cesaroni. & Shinnaga 



core JCMT 18354-0649S are detected (IWu et al 



201 



20061 : 



Klaassen. & Wilson 



2003 



2OO7I : 



Wu et al 



). Both infall and outflow motions in the massive 



observations (Liu et al 



20051) . and further confirmed by higher resolution 



20111). Although accretion-disk-outflow systems are found in high-mass 



star forming regions, there may be differences between low- and high-mass formation. 

The infall motion can be detected via "blue profile", a double-peaked profile with the 
blueshifted peak being stronger for optically thick lines and a single peak at the absorption part of 
optically thick lines for optically thin lines, whic h is caused by sel f absorption of the cooler outer 



infalling gas towards the warmer central region (|Zhou et al 



19931). In contrast, the "red profile" 
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where the redshifted peak of a double-peake d profile being 



suggested as indicators for outflow motions. 



s tronge r for optically thick lines is 



Mardones et al. 



(119971) defined the "blue excess" in a 



survey, E, as E = (Nb-N«)/N7- (Mardones et al. 1997), where Nj- is number of sources, and Nr 
mark the number of sources with blue and red profiles, respectively. The blue excess seems to be 
no significant differences among the low-mass co res in diff"erent evolutionary phases. However, 



using the IRAM 30 m telescope, Wu et al. 



(|2007i) found that UC Hii regions show a higher blue 
excess than their precursors, indicating fundamental diff"erences between low- and high-mass-star 
forming conditions. The searches need to be expanded. 



Located at a distance of 5.7 kpc (|Hofner et al.lll994|) . G9. 62-1-0. 19 is a well studied high-mass 
star forming region containing a cluster of Hii regions, which are probably at diff"erent evolutionary 



stages. Mu 



from A-I) (IGarav et al 



tiwavelength VLA observations hav e identified nine radio continuum sources (denoted 



in diameter) (|Garay et al 



1993 



1993 



esti et al 



Testi et al 



20001). and components C-I are very compact (< 5' 



20001) . As revealed in NH3 (4,4), (5,5) and CH3CN 



(J=6-5), co mponent F is a hot molecular core (HMC) and hence likely the youngest source in 



the region (Cesaroni et al 



19941 : 



Hofner et al. 



199^ 



1996bl) . G9.62-h 0.19 E is a young m assive 



star surrounded by a very small UC Hii region and a dusty envelope (IHofner et al. 



G9.62+0. 



Testi et al. 



9 D a small cometary UC Hii region excited by a BO. 5 ZAMS star (|Hofner et al. 



1996b L while 



1996bl : 



20001) . Both G9. 62-^0.19 E and G9.62-h0.19 D seem to be at a more evolved stage than 
G9. 62-1-0. 19 F. Thus G9.62 complex is an ideal sample to examine massive star forming activities 
including outflow and infall motions. 

Maser emissions of NH3, H2O, OH, and CH3OH, as well as the strong therm al NH3 emissi ons 



were detected along a narrow region with projected length 20" and width< 2" (|Hofner et al 



1994h . A possible explanation for this alignment is compression of the molecular gas by shock 
front originating fro m an even more evolved Hii region to the west of the star-forming front 



(IHofner et al. 



1994|) . High-velocity molecular outflows also have been detected in this region. 
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and G9.62+0.19 F is believed to be the driving sourc e (IGibb. Wyrowski. & Mundy 



Hofher. Wiesemeyer. & Henning 



2001 



SuetaL 



2004 



20051). However, most of previous work was 



carried out at low frequencies, probing low excitation conditions. To exam the hot dust/gas 
environment and dynamical processes in this region, higher resolution studies at high frequencies 
are needed. In this paper we report the results of the Submillimeter Array (SMAlj) observations 



toward G9.62+0.19 region at 860 /urn. 



2. Observations 

The observations of 09. 62+0. 19 with the SMA were carried out on July 9th, 2005 with 
seven antennas in its compact configuration at 343 GHz for the lower sideband (LSB) and 
353 GHz for the upper sideband (USB). The T^y, ranges from 210 to 990 K with a typical 
value of 380 K at both sidebands during the observations. The observations had two fields 
for the G9.62+0.19 complex to cover the entire region with emissions. One phase reference 
center was RA(J2000) = 18''06'^14.2P and DEC(J2000) = -20°31'46.2", and the other was 
RA(J2000) = 18'^06'"15.00'^ and DEC(J2000) = -20°31'34.20". Uranus and Neptune were 
observed for antenna-based bandpass calibration. QSOs 1743-038 and 1911-201 were employed 
for antenna-based gain correction. Neptune was used for flux-density calibration. The frequency 
spacing across the spectra band was 0.8125 MHz, corresponding to a velocity resolution of ~0.7 
km s~^. 



MIRIAD was employed for calibration and imaging dSault et al 



19951). The imaging was 



done to each field separately and the mosaic continuum map was made using a linear mosaicing 



' Submillimeter Array is a joint project between the Smithsonian Astrophy sical Observatory and 
the Academia Sinica Listitute of Astronomy and Astrophysics and is funded by the Smithsonian 
Institution and the Academia Sinica. 
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algorithm (task "linmos" in MIRIAD). The 860 yum continuum data was acquired by averaging 
over all the line-free channels in both sidebands. The spectral cubes were constructed using 
the continuum- subtracted spectral channels smoothed into a velocity resolution of 1 km s ^ 
Additional self-calibration with models of the clean components from previous imaging process 
was performed on the continuum data in order to remove residual errors due to phase and 
amplitude problems, and the gain solutions obtained from the continuum data were applied to the 
line data. The synthesized beam size of the continuum emission with robust weighting of 0.5 is 
2.76" X 1.88" (P.A.=21.4°). 



3. Results 



3.1. Continuum emission 



The 860 /im continuum image combining the visibility data from both sidebands is shown 
in Figure 1. Thre e sub-mm cores are detected. The known cm and mm continuum components 



(Testi et al. 



20001) of B. C. D. E, F, G, H, and I are marked by plus signs. Wate r masers 



(Hefner. & Churchwell 



1996ah are mark ed by open squares and methanol masers (INorris et al. 



1993|) by triangles. The near-IR sources (iPersi et al.ii2003i : ITesti et al, 



1998 



Linz et al. 



2005h are 



marked by filled circles. IRAC sources are taken from the database of Galactic Legacy Infrared 
Mid-Plane Survey Extraordinaire (GLIMPSE) □ and labeled with asterisks. The northern core is 
located at south-east of G9. 62-1-0. 19 C, and the middle core is associated with G9. 62-1-0. 19 E. 
The 860 /urn continuum emission at the southern core is concentrated on the hot molecular core 
G9.62+0.19 F and extends to G9.62+0.19 D in the south and to G9.62+0.19 G in the north. 



Gaussian fits were made to the continuum. The northern core seems to be a point-like source. 



%ttp://irsa.ipac.caltech.edu/data/SPITZER/GLIMPSE/ 
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The middle core is very compact with a deconvolved size of ~ 1.4". The southern core is found 
to be elongated from north to south with an average size of 2.4", containing at least three sources, 
D, F, and G. F is at its peak position. The peak positions, sizes, peak intensities and total fluxes 
of these three sub-mm cores are listed in Column 2-5 in Table 1 . The physical properties of these 
cores will be further discussed in section 4.2. 



3.2. Line emission 

Tens of molecular transitions including hot molecular lines CH3OH, HCOOCH3, and 
CH3OCH3 are detected t oward both the middle and southern sub-mm cores, indicating these two 



cores are hot and dense (|Oin et al 



120101) . Figure 2 presents the full LSB and USB spectra in the 
UV domain over the shortest baseline. The strongest lines are identified and labeled on the plots. 
Only HCN (4-3) and CS (7-6) line emissions are detected towards the northern core with our 
sensitivity. Thus we mainly focus on the middle and southern cores in this paper. The systemic 
velocities of 2.1 km s"^ for the middle core and 5.2 km s"^ for the southern core are obtained 
by averaging the V/^,- of multiple singly peaked lines. Six transitions of the thioformaldehyde 
(H2CS) and molecular transitions SO (87-77), CS (7-6), HC^^N (4-3) and HCN (4-3) are analyzed 
here, while the others will be discussed in another paper. We have made gaussian fits to the 
beam-averaged spectra, and present the observed parameters of these lines in Table 2. 



3.2.1. Line emission at the middle sub-mm core 

The integrated intensity maps of four transitions of H2CS towards the middle core are shown 
in the upper panels of Figure 3. From (a) to (d), the upper level energy of H2CS transitions varies 
from ~ 90 K to ~ 400 K. The H2CS emission is spatially coincident with continuum emission 
of the middle core very well. The Position- Velocity (P-V) diagram and first moment map of 
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HiCS (102,8-92,7) emission are presented in Figure 4. The P-V diagram is constructed across 
the peak of the continuum along the N-S direction. From P-V diagram two emission peaks are 
clearly revealed. The velocities of the two emission peaks are at 1 and 3 km s"^ with 1.5" spatial 
separation, indicating a velocity gradient in N-S direction. The first moment map also shows 
velocity changes in N-S direction. The small velocity gradient detected in H2CS (102,8-92,?) 
emission may indicate a disk with a low inclination along the line of sight, which requires further 
confirmation with higher angular resolution observations and other molecular line tracers. 

The spectra and integrated intensity maps of HC^^N (4-3) and SO (Sy-V?) are presented in 
Figure 5. The two spectra seem to be symmetric, and their cores are associated with that of the 
continuum emission very well. 

Figure 6 presents the spectra and P-V diagrams of HCN (4-3) and CS (7-6) emissions of the 
middle core. HCN (4-3) and CS (7-6) show asymmetric profile. The blue and red emission peaks 
of HCN (4-3) are around km s~^ and 6 km s~\ respectively. The blueshifted emission of CS 
(7-6) peaks around 1 km s~\ while the redshifted around 4 km s~^ We can see the blueshifted 
emission of both HCN (4-3) and CS (3-2) is always stronger than the redshifted emission and the 
absorption is also redshifted, which are blue profiles (see Sect. 1). Besides the "blue profile", 
some weak absorption dips are found around 10 km s"' in both the spectra and P-V diagrams of 
HCN (4-3) and CS (7-6), and further observations are needed to determine the properties of these 
absorption dips. In this paper we only pay attention to the "blue-profile" found in CS (7-6) and 
HCN (4-3) emission. 

The integrated intensity maps of HCN (4-3) and CS (7-6) towards the middle core are 
presented in Figure 7. The HCN (4-3) and CS (7-6) are associated with the dust emission. 
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3.2.2. Line emission at the southern core 

The integrated intensity maps of four transitions of H2CS at the southern core are shown in 
the lower panels of Figure 3. The upper level energy of H2CS transitions varies from ~ 90 K to 
~ 400 K from panel (e) to panel (h). As the upper level energy increases, the emission peak of the 
different transitions of H2CS moves from S-E to N-W, indicating a temperature gradient in the 
southern core. 

Averaged spectra of SO (Sy-Vy), HC^^N (4-3), HCN (4-3) and CS (7-6) at the southern core 
are presented in Figure 8. The spectra of SO (Sy-V?) and HC^^N are averaged over a region of 
4", while HCN (4-3) and CS (7-6) are averaged over a region of 6". SO (Sy-V?) emission has a 
total velocity extent of larger than 20 km s~' . From gaussian fit to the spectrum, the peak velocity 
of SO emission is 5.1 + 0.1 km s"\ coincident very well with the systemic velocity 5.2 km s"'. 
HC^^N (4-3) has a velocity extent of about 15 km s~^ The velocity extents of CS (7-6) and HCN 
(4-3) are as high as 40 km s~' and 60 km s~\ respectively. Emission wings are clearly detected 
from the spectra of the four lines. A "red-profile" is significantly exhibited in the spectra of CS 
(7-6) and HCN (4-3), of which the redshifted emission is always stronger than the blueshifted 
emission with an absorption dip at the blueshifted side of the systemic velocity (5.2 km s"'). This 
profile is caused by absorption of the colder blueshifted gas in front of the hot core, indicating 
outflow motion s. The "red-profile" is consistent with tha t detected using single dish observations 



(see Figure 6 of 



Hofner. Wiesemeyer. & Henning 



mom . 



The integrated intensity maps of HC^^N (4-3) and SO (87-77) at the southern core are 
presented in Figure 9. To avoid the influence of outflow motions, both the maps are integrated 
from 2 km s"' to 8 km s~'. Both the emission of HC^^N (4-3) and SO (87-77) coincides with the 
cm/mm component F, and extends from D to G. 

As shown in the left panels of Figure 10, the high velocity gas of HC^^N (4-3) and SO (87-7?) 
can be identified by the vertically dashed lines in the P-V diagrams. For SO (87-77), we integrate 
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from -4 km < V < km for the blue wing and 10 km < V < 14 km for the red wing, 
and present the contour map in (c) of Figure 10. For HC^^N (4-3), only the red wing emission is 
presented in (d) of Figure 10. The high velocity emission of both HC^^N (4-3) and SO (Sy-Vy) 
is associated with core F, indicating that core F is the driven source of the outflow. The blue and 
red wings of SO (S?-??) overlap to a large extent in the contour maps, and hence the molecular 
outflow revealed by SO (Sy-Vy) is observed close to its flow axis. 

Figure 1 1 presents the channel maps of CS (7-6) emission. The redshifted high- velocity gas 
seems to be elongated from north-east to south-west, while the blueshifted high- velocity gas from 
north to south. The high velocity gas revealed by CS (7-6) is also very obvious in the P-V diagram 
in Figure 13(d). As shown in the P-V diagram, the blueshifted high- velocity gas extends about 8" 
from north to south. The high-velocity emission integrated over the wings (-12 km s~^ < V < -5 
km s~^ for the blue wing and 15 km s~' < V < 22 km s"' for the red wing) is presented in Figure 
13(e). 

Figure 12 is the channel maps of HCN (4-3) emission. The maximum of absorptions appears 
at around km s~^ The redshifted high-velocity gas seems to be elongated from west to east, 
while the blueshifted high-velocity gas from north to south. At very high velocity channels (V< 
-16 km s"^), the blueshifted emission is totally located at south-east. By comparing the channel 
maps and P-V diagrams (see Figure 13) of HCN (4-3) and CS (7-6) at velocity intervals -12 
km s~^ < V < -5 km s~^ and 15 km s~^ < V < 22 km s~\ we find similar structures in CS (7-6) 
and HCN (4-3) emissions. The high-velocity emission of HCN (4-3) integrated from -12 km s~^ 
to -5 km s"^ for the blue wing and from 15 km s"^ to 22 km s"^ for the red wing is presented in 
the panel (b) of Figure 13. As of CS (7-6), the blueshifted gas revealed by HCN (4-3) is elongated 
from north to south with the emission center located between G9. 62+0. 19 F and G9.62 -1-0.19 D, 
while the redshifted gas from north-east to south-west. Two clumps are found in the blueshifted 
high- velocity emission of HCN (4-3), which locate at north-west and south-east of F, respectively. 
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In order to reveal the very high velocity emission traced by HCN (4-3) but not CS (7-6), we 
integrate over the wings at much higher velocities (-20 km s"' < V < -13 km s"' for the blue 
wing and 23 km s"^ < V < 39 km s~^ for the red wing), and present the integrated emission map 
in Figure 13(c). The redshifted emission is elongated from north-east to south-west with the 
emission center located between G9. 62-1-0. 19 F and G9. 62-1-0. 19 G, while the blueshifted emission 
center located between G9. 62-1-0. 19 F and G9.62-l-0.19 D. It is clearly seen that the high velocity 
gas traced by SO (87-77), CS (7-6) and HCN (4-3) have different spatial distributions, which 
should be caused by the complicated interactions between the outflow and the ambient gas. It 
may also indi cate a change o f the outflow axis. The ch ange of outflow axis is also found in IRAS 



20126+4104 (Suetal 



20071) and JCMT 18354-0649S ( Liu et al. 



2011 ). 



From the integrated emission maps of SO (87-77), HCN (4-3) and CS (7-6) high velocity gas, 
it is clearly seen that G9. 62-1-0. 19 F is located at the middle of the redshifted and blueshifted lobes, 
suggesting G9. 62-1-0. 19 F is the outflow driving source. 



4. Discussion 



4.1. Rotational temperature of H2CS transitions 



Six transitions of H2CS have been detected in the middle and southern cores, enabling us to 
estimate the rotational temperature. Under the assumptions that the gas is optically thin under 
local thermodynamic equilibrium and the gas emission fills the beam, the rotation temperature 
and beam- averaged column density can be estimated using the Rotational Temper ature Diagram 



(RTD) by (|Cummins. Linke. &Thaddeus 



19861 : 



Turner et al. 



19911: 



Liu, et al. 



2OO2I) 



In(^) = In(^) ln[2.04 x 10^° 



J I(Jy beam ^)dv(km s ') 



-] 



Qro/ Trot ea9tiarcsec^)gigKV^(GHz^)Siu\debye^) 
where N„ is the observed column density of the upper energy level, gi, is the degeneracy factor 
in the upper energy level, N7- is the total beam-averaged column density, Q^f is the rotational 



(1) 
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partition function, E„ is the upper level energy in K, T^ot is the rotation temperature, J I dv is 
the integrated intensity of the specific transition, 0„ and Oh are the FWHM beam size, gK is the 
K-ladder degeneracy, g/ is the degeneracy due to nuclear spin, v is the rest frequency, and S is 
line strength and the permanent dipole moment. For H2CS, the interchang eable nuclei are sp in 



^, leading to ortho - and para- forms with g/ equaling | and 7 , respectively (Blake et al 



Turner et al. 



199 ih . The partition function Qrot of H2CS is (Blake et al 



1987 ) 



19871 : 



Qrot = 2[ 



n{kTrotf 



(2) 



h^ABC 

where k and h are the Boltzmann and Planck constants, respectively, and A, B, and C are 
the rotation constants. Thus the rotation temperature T^ot and total column density can be 
estimated by least-squares fitting to the multiple transitions. We applied the RTD method towards 
D, E, F, G (see Figure 14), and the fitting results are listed in the second and third columns of 
Table 3. The rotational temperature of the middle core (E) is 83+21 K. In the southern core, the 
rotational temperature estimated decreases from G (91 K) to F (83 K) and D (43 K), suggesting 
the temperature gradient in the southern core. The total column density of H2CS ranges from 
1.3x10^5 (G) to 3.8x1015 cm-2 (D). 

However, the filling factor and the optical depth correction were not taken account of in 
the RTD method. To investigate their effect vy e applied the Population Diagram (PD) analysis 



(IGoldsmith. & Langer 



19991 : 



Wang et al.l 



2010f) . In the PD analysis, we have 

) 



In(^) = In(^) - ^ + ln(/) - ln( 



El' 

Trot 



Qrot Trot 1 £ 

where A^„ is the inferred column density of the upper energy level from the PD analysis, f is 
the source filling fac tor and t is the optical depth. The optical depth r can be expressed by 



(3) 



(IRemijan et al. 



2004 ) 



%TfSl/v Nt .Ik. 

T = e ^rol 

rot 



(4) 
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where Av is the FWHM line width. Under LTE, the upper-level populations, Nu, can be 
predicted according to the right-hand side of Equation (3) for a given set of total column density, 
Nj-, rotational temperature. Trot, and source filling factor, f. The expected A^^ were evaluated for 
the parameter space of Trot = 10-500 K, Nr = lO^^^-lO^^ cm-^, and f between 0.01 and 1.0. To 



where 5 Nu is the 1 cr error of observed upper-state column density. Although the is a 
good representation of the goodness of fit, the parameter set with the lowest rn^Y ^ot actually 
represent physical parameters very well due to the uncertainties of the observed data. In order to 
find a representative parameter set, we compute a weighted mean and standard deviation for all the 
parameters, with the weights being the inverse of the;^^. All the parameter sets where the inferred 
upper- level population corresponds with the observed upper-level population Nu within 3 cr are 
used to compute the weighted means and standard deviations. The derived rotational temperature, 
total column density and filling factor of each component are list in the [3-5] columns of Table 3. 
The inferred optical depths of each line transition are listed in the last six columns of Table 3. The 
rotational temperatures of D, E, F, G are estimated to be 42+34, 92+74, 51+23 and 105+37 K, 
respectively. A temperature gradient in the southern core is also revealed as in the RTD method. 
The four components D, E, F, G has similar total column densities as high as 4x10^^ cm~^, about 
an order of magnitude higher than those obtained from RTD method, which are mainly due to the 
small source filling factor (<0.5). The optical depths of H2CS (10o,io-9o,9) at the four components 
are all much larger than one, while the other transitions are always optically thin except H2CS 



compare the observed Nu and the inferred Nu, we calculate thex^ as: 




(5) 



(102,9-92,8) line at G. 
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4.2. Core properties 



In the optically thin case, the total dust and gas masses of the three sub-mm cores can be 



19831) . where Sy is the flux at 860 



obtained with the formula M = S yD^ / KyRBy(Td) ([Hildebrand 
fim, D is the distance, R=0.01 is the mass ratio of dust to gas, and Ky is dust opacity per unit dust 
mass. By{Td) is the Planck function at a dust temperature of T^. We assume that equals the 
rotational temperature of H2CS. For the northern core, since only CS (7-6) (upper energy = 65.8 



K) and HCN (4-3) (E„ = 42.5 K) exhibit strong emission lines, we assume 



SuetaL 



to be 50 K. Together 



(l2005h extrapolated 



with the measurements at centimeter and millimeter wavelengths, 
the ionized gas emission at mm/submm wavelengths, and found that the 0.85 mm continuum 
associated with components D, E, and F are dominated by thermal dust emission. They have 
derived opacity index of components E and F to be 1.2, and 0.8, respectively. For the northern 
sub-mm core, /3 = 1.5 is assumed. Using the above dust opacity indexes, we adopt Ky=2.Q, 1.8, 
and 1.5 cm^g"^ for the northern, middle and southern cores, respectively (jOssenkopf & Henning 



1994h . At the distance of 5.7 kpc , we get the total dust and gas masses for these three cores, and 
list all the parameters in Table 1 . The deduced masses for the northern, middle and southern cores 
are 13, 30, 165 Mq, respectively. The column density of H2 are 1.2x10^"* and 2.1x10^"* cm~^ for 
the middle and southern sub-mm cores, respectively. 



4.3. Infall properties in the middle core 



In the middle core, both CS(7-6) and HCN(4-3) emission exhibi ts "blue profile" featur e. 



indicating infall motions of the gas envelope 


towarc 


the central star ( 


Keto. Ho.& 


Hasc 


lick 




1988; 


Zhou et al. 


1993; 


Zhan 


1. Ho. & Ohashi 


1998; 


Wu & Evans 


2OO3I 


Wu et al. 


2005, 


2007; 


Fuller. Williams. & Sridharan 


2005; 


Wvrowski 


2007 


: Sun. & Gao 


20081). 


The velocity difference 



(0.9 km s ^) between the absorption dip in CS (7-6) spectrum (3 km s ^) and the systemic velocity 
(2.1 km s"') is taken as the infall velocity y„,. Since both HCN (4-3) and CS (7-6) emissions are 
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not resolved towards the middle core, we simply take the dust core size as the radius of the infall 
region, which may underestimate the infall rate derived below. The kinematic mass infall rate 
can be calculated using dM/dt=4nR^^nmVin. n=1.5xlO'cm"^ is the number density of this dust 
core. Taking Helium into account, the mean molecular mass m is 1.36 times of H2 molecule mass. 
The infall rate calculated is 4.3 x 10"^ Mo-yr"^ For comparison, the Vi„ from pure free-infall 
assumption is also derived with the formula Vf^ = 2GM/Ri„. The pure free-infall velocity is 
Vin = 3.6 km s"' and thus the "gravitational" mass infall rate is 1.7 x 10"^ Mo-yr"\ which is larger 
than the kinematic infall rate. 



4.4. Outflow properties in the southern core 



4.4.1. Shock chemistry in the outflow region of the southern core 



different outflow regions (Bachiller et al. 



2004 



Codella et al. 



1997 



Observations have suggested t hat there ar e important diffe renc es in molecular abundances in 

:Choietal.l2004: 



Jorgensen. Schoier. & van Dishoeck 



20051). Si gnificant abundance enhancements are found in the shocked regi on 



for sulfur-bearing molecules (|Bachiller et al 



1997 



Jorgensen. Schoier. & van Dishoeck 



and the abundance of HCN in outflow regions is related to atomic carbon abundance (IChoi 



2004) 



20021) 



However, previous studies of the chemical impact of outflows are confined to the well coUimated 
outflows around C lass sources, while such studies especially high resolution studies on mass ive 



outflows are rare (iBachiller et al. 
2007). 



19971 : 



Jorgensen. Schoier. & van Dishoeck 



2004 



Arce et al 



A red and bright IRAC source is found to be associated with the southern core. The 
magnitudes of the IRAC source at 3.6 //m, 4.5 yum and 5.8 yum are 10.102 + 0.093, 8.361 + 0.108 
and 7.778 + 0.302 mag, respec tively. The [3.6-4.5 1 color is as large as 1.74, indicating shocked 



emission in the southern core (ITakami et al 



2010|) . Maser emissions of NH3, H2O, OH, and 
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CH3 OH, as well as the 



gas (|Hofner et al. 



strong thermal NH3 emissions also uncover the existence of the shocked 



1994|) . Outflows can be reveale d from shocked emission probe d by the 



2008 



Takami et al. 



2010h . Thus the 



strong and extended emission at the 4.5 pim band (lOiu et al.l 
massive outflow in the southern core of G9.62 complex provides an ideal sample to study shock 
chemistry. 

The fractional abundance of a certain molecule is defined as;^f = Nt/Nhj, where Nt is the 
total column density of a specific molecule and is the H2 column density. Assuming that the 
gas is optically thin and the emission fills the beam, the beam-averaged total column density of a 
specific molecule can be obtained from: 



Nt = 2.04 X 10 



20 



J I{Jy beam ^)dv(km s ^)Qrote^"'^" 



(6) 



6a9},{arcsec'^)gigKV^{GHz^)S ^^debye^) 

Assuming that T,.o, of HC'^N equals to that of H2CS and the gas is optically thin, Nt of 
HC^^N is calculated to be 3.0 x 10^^ cm~^ at the core region. At the galactocentric distance of 3 kpc 
for G9.62+0.19 (Scoville et al.l987,Hofner et al.l994), the abundance ratio [^^NJA^^N] « 350 
(Wilson and Rood. 1994). Thus the total column density of HCN at the core region should 
be 1.1 X 10'^ cm^^. Therefore, the fractional abundance of HCN relative to H2 at the core 



region is 5.2 x 



(Bachiller et al. 



^. HCN appears to be greatly enhanced in the outflow regions of the LI 157 



19971), wh ile has similar abundances in the outflow re gion and the ambient cloud 



of NGC 1333CIRAS 2A (Ijorgensen. Schoier. & van Dishoeck 



2004 . Owing to the lack of a 



direct estimation of the H2 column density towards the outflow region, the fractional abundance 
of HCN in the outflow region is also assigned to 5.2 x 10"^ in calculating the outflow parameters. 
Since the HC'^N emission traces outflowing gas at much lower velocity than HCN, perhaps HCN 
could be more enhanced in the high velocity component. With the possibility of higher opacity 
and the lack of direct H2 column density measurement, the derived fractional abundance perhaps 



is a lower limit anyway. 



Su et al. 



(20071) estimate an HCN abundance of ~ 1 - 2 x 10 Mn the 
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massive outflow lobes of IRAS 20126+4104, which is comparable to our estimation here. 

Since the blueshifted outflow gas traced by CS (7-6) and HCN (4-3) suffers self-absorption, 
the abundance ratios among SO (8? - 7?), CS (7-6), and HCN (4-3) were inferred from the 
beam-averaged spectra taken from the redshifted outflow lobe. The abundance ratio as a function 
of flow velocity (the outflow velocity relative to the systemic velocity) of [CS/SO] is obtained 
assuming five different excitation temperatures in the left panel of Figure 15. It can be seen 
that the abundance ratio of [CS/SO] increases with the excited temperature. At each excitation 
temperature, the abundance ratio of [CS/SO] has lower values at flow velocities less than 6 km s~\ 
and higher values when W flow larger than 8 km s"\ whereas the abundance ratio seems to 
be constant at flow velocities between 6 km s ' and 8 km s~'. There are two reasons for the 
lower abundance ratio when N flow < 6 km s first, the flux missing of CS (7-6) due to the 
interferometer is more serious than SO (8? - 7?); second, CS (7-6) may be more optically thick at 
lower flow velocities than SO (87 - 7?). As shown in the P-V diagrams, the emission region of 
CS (7-6) is much larger than SO (87 - 7?) at high velocities. The higher abundance ratio when 
V//OH, > 8 km s ' is due to the smaller filling factor of SO (87 - 77) emission. We propose the 
mean observed value between 6 km s^ and 8 km s^ can represen t the actual abu ndance ratio of 



[CS/SO]. Assuming a typical excitation temperature of T^v =30 K dWu et al 



Nilsson et all (2000h find that the 



2004), the abundance 

T 



ratio of [CS/SO] at the redshifted lobe is inferred as 0.7. 
[SO/CS] abundance ratios are strongly enhanced in the Orion A and NGC 2071 outflows where 
the [SO/CS] ratios are estimated to be about 24 and 2.2, respectively. However, the [SO/CS] 
abundance ratio in the outflow of G9. 62-1-0. 19 is found to be 1.4, much lower than that found in 
Orion A outflow. 

As shown in the right panel of Figure 15, the abundance ratio of [CS/HCN] decreases linearly 
with the flow velocity. To avoid the missing flux difficulty, the abundance ratio is calculated at 
high flow velocities larger than 7 km s~^ The decreasing of the abundance ratio with velocity is 
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because that the emission region traced by CS (7-6) is always smaller than HCN (4-3), leading to 
smaller filling factor for CS (7-6), which can be verified easily by comparing the channel maps 
between CS (7-6) in Figure 11 and HCN (4-3) in Figure 12 at high velocities. We fitted the 
observed data with a linear function, and adopted the value at flow velocity of 10 km s ' as the 
actual abundance ratio of [CS/HCN] in the outflow region, which is [CS/HCN] = 1.2. Since HCN 
fractional abundance is 5.2 x 10"^, the fractional abundances of CS and SO are deduced to be 
6.2 X 10"'' and 8.9 x lO"'', respectively. 



4.4.2. Properties of the bipolar-outflow traced by SO (87 - 7? j emission 



The SO (87 - 7?) emission in the southern core shows line wings, suggesting outflow 
motions. From the integrated intensity map in Figure 10(c), we find the outflow lobes revealed 
by SO (87 - 77) emission peak at diff erent position with different positio n angle compared with 



previously reported H7S (2? - 2i 1 ) (IGibb. Wvrowski. & Mundy 



(|Hofner. Wiesemeyer. & Henning 



20041) and HCO+ (1-0) data 



20011). But in the same sense, the blue- and red-lobes revealed 
by SO overlap to a large extent as well as HCO^ ( 1-0) and H7S (l^o - 2i 1) data, consisten t with 



the argument of the outflow being viewed pole-on (|Hofner. Wiesemeyer. & Henning 



200 ih . 



The total mass of each outflow lobe is given by: 



M 



flow 



1.04 X IQ-^D^ 



Sydv 



(7) 



where M//o»., D, Sy, and r are the outflow gas mass in M©, source distance in kpc, line flux 
density in Jy, relative abundance to H2, and optical depth. The other parameters have the same 
units as in equation (1). The fractional abundance of SO is taken as 8.9 x 10"^ (see Sec. 4. 4.1). 
Assuming an excitation temperature of 30 K and the outflowing gas is optically thin, the inferred 
outflow masses are 13 M© for each of red and blueshifted lobes. Thus, the momentum can be 
calculated hy P = 2M(v)dv, and the energy hy E = Yj ^M(v)v^dv, where v is the flow velocity. 
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The derived parameters are listed in Table 4. The momentum and energy of the red lobe are 82 
Mg-km s"^ and 5.4 x 10"^^ erg. For the blue lobe, the momentum and energy are calculated to be 
86 Mg-km s~' and 5.8 x 10"^^ erg. The dynamical timescale tj„i is estimated as RfVchar, where 
R (~ 0.06 pc) is adopted as the mean size of the outflow lobes assuming a coUimation factor of 
unity, and Ychar (~ 5.5 km s"') is assumed as the mass weighted mean velocity. Thus, the dynamic 
timescale is estimated to be 1 x 10"* year, which may be underestimated due to the uncertainty 
of the outflow scale. The mechanical luminosity L, and the mass-loss rate M ar e calculated as 



L=E/ t, M = PI{tVw), where the wind velocity V», is assumed to be 500 km s ^ ([Lamers et al. 



19951). The mechanical luminosity L and the total mass-loss rate are estimated to be 9.3 Lq and 
3.6 x 10"^ Mo-yr"^ respectively. 



4.4.3. Very high-velocity gas detected in CS (7-6) emission 

The CS (7-6) emission at the southern core shows "red-profile" with wide wings. We take 
6.2 X 10"^ as the fractional abundance of CS relative to H2 along the outflow lobes. Assuming 
Tex = 30 K, we derive the parameters for the CS outflow (Table 4) with the same method used for 
SO (87 - Vv). The outflow masses at very high velocities {\fiow > 10 km s"^) are 3.7 Mq and 
5.5 Mq for the blueshifted and redshifted lobes, respectively. The momentum and energy of the 
blueshifted lobe at very high velocities are calculated to be 47 Mg-km s"' and 6.0 x 10"*^ erg. For 
the redshifted lobe, the momentum and energy at extremely high velocities are calculated to be 68 
Mo-km s~^ and 8.7 x lO'^^ erg, which are similar to the blueshifted lobe. 



4.4.4. Very high-velocity gas detected in HCN (4-3) emission 

As discussed before, HCN (4-3) has a velocity extent of at least 60 km s"\ which traces 
extremely high- velocity (EHV) gas. Adopting an excited temperature of 30 K, and an HCN-to-Ha 
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abundance ratio of 5.2 x 10"^, the parameters of the outflow are calculated and listed in Table 
4. The outflow mass at very high velocities {Vfiow > 10 km s"') are 5.2 Mq and 17.6 Mq for 
the blueshifted and redshifted lobes, respectively. The momentum and energy of the blueshifted 
lobe at very high velocities are 85 Mg-km s"' and 1.4 x 10"^^ erg. For the redshifted lobe, the 
momentum and energy at very high velocities are 294 Mo-km s"' and 5.5 x 10"^^ erg, which are 
larger than the blueshifted lobe. 



4.4.5. Mass-Velocity diagrams 



A broken power law, dM(v)/dv cc y y usual 



2004 



stellar objects (IChandler et al. 



Oiu et al. 



2007 



19961 : 



Lada. & Fich 



V exhibits in molecular outflows near young 



19961 : 



Ridge. & Moorell2001 



Su. Zhang. & Lim 



20091) . The slope, y, typically ranging from 1 to 3 at low outflow velocities, 
and often steepen s at velocities larger than 10 km s^' — with y as large as 10 in some cases 



(Arceetal 



20071) . Assuming optically thin, the mass-velocity diagrams of the outflow at the 
southern core of G9.62+0.19 complex are shown in Figure 16. SO (87 - 7?), CS (7-6), HCN (4-3) 
results were all used in the mass spectra. We calculate the outflow mass traced by CS (7-6) and 
HCN (4-3) from V//„„, of 10 km s"' to avoid the absorption of the spectra. Instead of broken 
power law appearance, the mass-velocity diagram of blueshifted lobe can be well fitted by a 
single power law with a power indexes of 2.28 + 0.23. The mass-velocity diagram of redshifted 
lobe can be well fitted by a single power law with a power indexes of 1.70 + 0.17 even though 
the mass drops more rapidly after 25 km s~'. As marked by the dashed ellipse in the right panel, 
the outflow mass revealed by CS (7-6) is much lower than that revealed by HCN (4-3) at very 
high velocities. Despite the CS data, the mass- velocity diagram of redshifted lobe at velocities 
smaller than 25 km s"' can be fitted by a single power law with a much smaller power indexes of 
1.08 + 0.09. However, no significant slope changes are found in both the red- and blue-shifted 
lobes of the outflow at the southern core, which are very diff'erent from those previous works. 
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4.5. Different evolutionary stages of the three dust cores 



The northern core has the smallest diameter and mass among the three cores. It seems 
likely to be a point source after deconvolution. It is located south of the nominal radio UC Hii 
region G9.62+0.19 C. In this region, eight near-IR sources are detec ted in a diffuse near-IR 



20031). The reddest one c7 



nebulosity at the west of the radio emission peak (iPersi et al.l 
(18''06"'14.34%-20°31'25.0") is located within 1" of the radio peak, while the faintest one c8 
(18''06™14.42',-20°3r27.4") seems to be associated with the sub-mm core detected in SMA 
observation. Source c8 is too faint to be detected even at H band and also shows no emission 
at 12.5 yum. In contrast to the bright, rich molecular spectrum forest in the middle and southern 
sub-mm cores, the northern sub-mm core lacks strong molecular emissions. There is also no other 
early star forming signature such as masers associated with it. Since it is with near-IR emission 
and at the edge of the UC Hii region G9. 62-1-0. 19 C, the northern core may be just a remnant core 
in the envelope of UC Hii region G9. 62-1-0. 19 C, which needs further observations. 



The middle core is asso ciated with the hyper-compact Hii region G9. 62-1-0. 19 E (iGaray et al. 



1993 



Kurtz. & Franco 



20021). OH. H^O , and NH^ (5,5) ma s ers have been detected near th e radio 



emission peak (Forster & Caswel 



19891 : 



Hofner et al. 



19941 : 



Hofner. & ChurchwelL 1996a). Peri- 



odic c 



20091 : 



ass II methanol masers are also found in G9. 62-1-0. 19 



Goedhart. Gavlard. & van der Walt 



2005 



Norris et al 



van der Walt. Goedhart. & Gavlard 



19931). Methanol masers are believed 



to be a good tracer of y oung massive star form ing regions at stages earlier than relatively 



evolved UC H 



E ( Persi et al. 



II reg ions (|Longmore et al 



2007 ). No infrared source coincides with G9. 62-1-0. 19 



20031). Hot molecular CH3CN lines are detected in this region, and a kinematic 



tempe rature of Tyt = 108 K was obtained from CH3CN emission with LVG model (jHofner et al 



1996b|) . which is coincident with the rotational temperature (T^r = 92 K) obtained from H2CS 
emission. A spectra forest including hot molecular lines, such as CH3OH, is detected towards 
G9. 62-1-0. 20 E, suggesting this core is in a hot phase. Infall motions are traced by CS (7-6) and 
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HCN (4-3) lines, indicating active star fonning in this region. All above suggest that G9. 62+0. 20 
E is forming a massive young star. 

The 860 //m dust emission of the southern core peaks at G9. 62+0. 19 F, and extends from 
north to south. A hump structure is found to the southeast of the emission peak, indicating 
another possible sub-mm core. The previously recognized mm/cm cores (G9. 62+0. 19 D, G) are 



at the edges of the southern core. G9. 62+0. 19 G is a weak radio source (|Testi et al 



G9. 62+ 0. 19 D is consistent with an isothermal UC Hii region excited by a BO. 5 star (|Hofner et al. 



2000I). while 



1996bh . Weaker radio emission was fou nd at core F. H7Q and 



whole sub-mm core from north to south (Forster & Caswell 



1989 



masers are found across the 



Hofner. & Churchwell 



A near-IR source with large NIR excess is found to be assoc iated with G9. 62+ 0. 19 F (|Testi et al 



1998 



Persi et al. 



Linz et al. 



1996a) 



(20051) found four near-IR 



20031) . With higher resolution observations 
objects (F1-F4) in this core. F4 is with little emission at K band but becomes redder at longer 
wavelengths, which seems to correspond to the bright IRAC source with large excess at 4.5 /im. 
This object is the dominating and closest associated source of core F. Core F is also confirmed to 
be the driving source of an active outflow. All of above imply that G9. 62+0. 19 F is a very young 



massive star forming region. 



4.6. Blue excess in high-mass star forming regions 



Wu et al. (2007) found that UC Hii regions s how a higher blue excess than UC Hii precursors 



with the IRAM 30 m telescope. 



Wvrowski et al, 



(l2006h also detected large blue excess in UC Hii 



regions. "Blue profile" was detected with CS (7-6) and HCN (4-3) lines in UC Hii region 
G9. 62+0. 19 E, while "red profile" in hot molecular core G9.62+0.19 F, which coincides with their 
argument. The detection of infall signature in G9. 62 +0.19 E also coincides th e interpretation that 



material is still accreted during the UC Hii phase (jWuetal 



2007; 



Keto 



20021). Around younger 



cores, the outflow is more active and cold than UC Hii regions, which leads to more "red profile". 
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While in UC Hii regions, the outflows become weak. The surrounding gas of UC Hii regions is 
thermalized and the temperature gradient towards the central star is more likely to cause "blue 
profile", which results in the higher blue excess than UC Hii precursors. 



5. Summary 

We have observed the G9. 62+0. 19 complex with the Submillimeter Array (SMA) both in the 
860 fum continuum and molecular lines emission. The main results of this study are as follows: 

1. Dust continuum at 860 fim reveals three sub-mm cores in G9. 62+0. 19 star forming 
complex. With H2CS as the rotational temperature prober, the temperatures of E and F are 
estimated to be 92+74 and 51+23 K, respectively. The mass calculated are 13, 30, and 165 Mq for 
the northern, middle and southern core. 

2. In the middle core, HCN (4-3) and CS (7-6) spectra exhibit infall signature. The infall rate 
calculated is 4.3 x 10"^ Mg-yr"'. The detection of infall signature in G9. 62+0. 19 E coincides the 



interpretation that material is still accreted after the onset of the UC Hii phase (|Wu et al. 



2007). 



3. In the southern core, high-velocity gas is detected in SO (85 - 7?), CS (7-6) and HCN (4-3) 
lines. A bipolar-outflow with a total mass about 26 Mq and a mass-loss rate of 3.6 x 10"^ Me-yr"' 
is revealed in SO (88 - 7?) line wing emission. G9.62+0.19 F is confirmed to be the driving source 
of the outflows in the southern sub-mm core. The abundance ratios of [CS/SO] and [CS/HCN] 
in the outflow region are found to be 0.7 and 1 .2, respectively. The abundance ratio [CS/HCN] 
decreases with the flow velocity, indicating smaller outflow regions revealed by CS (7-6) than 
that revealed by HCN (4-3). The mass-velocity diagrams of the blueshifted and redshifted 
outflow lobes can be well fitted by a single power law. The power indexes for the blueshifted and 
redshifted lobes are 2.28 + 0.23 and 1.70 + 0.17. No significant slope changes are found in the 
mass-velocity diagrams. 
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4. The evolutionary sequence of the cm/mm cores in this region are also analyzed. The 
northern core may be just a remnant core in the envelope of UC Hn region G9.62+0.19 C, which 
needs further observations. The middle core (G9. 62+0. 19 E) is in a hyper-compact Hii region. 
Core G9. 62+0. 19 F is confirmed to be a hot molecular core. 

5. The detection of blue profiles at the hyper-compact Hii region E and the red profiles at the 
hot molecular core F supports the results of single-dish observations that UC Hn regions have a 
higher blue excess than their precursors. 
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Fig. 1. — The 860 /urn continuum emission image. The contour levels are from 0.03 Jy beam~^ 



(3cr) in steps of 0.06 Jy beam (6cr). The known cm and mm continuum c omponents (iTesti et al. 



2000h of B, C, D, E, F, G, H, and I are marked by plus signs. Water masers (IHofner. & Churchwell 



1996ah are mark ed by open squa r es and methanol masers (iNorris et al. 



near-IR sources (Persi et al..2003 



Testi et al 



19981 : 



Linz et al 



1993h by triangles. The 



2005|) are marked by filled circles. 



IRAC sources are marked with asterisks. 
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Fig. 2. — The full LSB and USB spectra in the UV domain over the shortest baseline. The strongest 
lines are identified and labeled on the plots. 
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RA (J3000) 

Fig. 3. — Integrated intensity maps of four transitions of H2CS at the middle (upper panels) and 
southern cores (lower panels). The known cm and mm continuum components are marked by plus 
signs as in the continuum map. The contour levels in all the panels are from 3cr in steps of 3cr. The 
rms levels are 0.3, 0.3, 0.3 and 0.2 Jy beam^^-km s"^ for H2CS (10ojo-9o,9) in panels (a) and (e), 
H2CS (102,8-92,7) in panels (b) and (f), H2CS (103,7-93,6) in panels (c) and (g), and H2CS (IO5-95) 
in panels (d) and (h), respectively. 
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Fig. 4. — The P-V diagram (left) and First moment map (right) of H2CS (102,8-92,?) emission at 
the middle core, (a) The contours of the P-V diagram are from 0.6 to 1.4 in steps of 0.2 Jy beam"^ 
(Icr). (b) contour plot of H2CS (102,8-92,?) integrated intensity image overlayed on the first moment 
map. The contours are from 0.9 (3cr) in steps of 0.9 Jy beam"^-km s~^ The First moment map is 
constructed from the data after imposing a cutoff of 3cr. 
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Fig. 5. — Spectra and integrated intensity maps of HC'^N (4-3) (upper panels) and SO (87 - 7?) 
(lower panels) at the middle core. The systemic velocity is marked with the thick vertical dashed 
lines at the spectra panels. The known cm and mm continuum components of C, and E are marked 
by plus signs at the integrated maps as the continuum map. (a) the beam-averaged spectrum of 
HC^^N (4-3) at E, (b) the integrated intensity map of HC^^N (4-3). The contour levels are -1.2 
(6cr), 1.2, 2.4, 4.2, 6.6, 9.6 Jy beam^^-km s~\ (c) the beam-averaged spectrum of SO (87 - Ij) at 
E. (d) the integrated intensity map of SO (87 - 7?). The contour levels are -1.2 (6cr), 1.2, 2.4, 4.2, 
6.6, 9.6, 13.2, 17.4, 22.2 Jy beam^^-km s'K 



-35- 



CO 
CD 



(0 

c 

CD 




5 10 15 
Velocity (km/s) 



o 

CD 

to 
o 



CD 
CO 





5 10 15 
Velocity (km/s) 



20 



Fig. 6. — Beam-averaged spectra and Position- Velocity (P-V) diagrams of HCN (4-3) (upper pan- 
els) and CS (7-6) (lower panels) at the middle core. The P-V diagrams are cut along a position 
angle of 0°. (a) the beam-averaged spectrum of HCN (4-3) at E, (b) the P-V diagram of HCN 
(4-3). The contour levels are -1.5 (Scr), -0.9, 0.9, 1.5, 2.1, 2.7, 3,3, 3.9, 4.5 Jy beam^^ (c) the 
beam-averaged spectrum of CS (7-6) at E. (d) the P-V diagram of CS (7-6). The contour levels are 
-1.5 (Scr), -0.9, 0.9, 1.5, 2.1, 2.7, 3,3, 3.9, 4.5 Jy beam^^ 
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Fig. 7. — Integrated intensity maps of HCN (4-3) (left) and CS (7-6) (right) at the middle core. 
The contour levels in both maps are from 1.5 (5cr) in steps of 3 Jy beam"' -km s~'. HCN (4-3) is 
integrated from -3 to 7 km s"\ while CS (7-6) from -1 to 6 km s"' 



-37- 




-10 10 20 -20 20 40 

Velocity (km/s) 



Fig. 8.— Averaged spectra of SO (87 - li) (upper-left), HC^^N (4-3) (lower-left), HCN (4-3) 
(upper-right) and CS (7-6) (lower-right) at the southern core. The spectra of SO (87 - 77) and 
HC^^N (4-3) are averaged over a region of 4", while HCN (4-3) and CS (7-6) are averaged over a 
region of 6". HCN v2 = I (4-3) emission is marked by the arrow in the upper-right panel, which 
can be clearly distinguished from the red wing of HCN (4-3). 
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Fig. 9. — The integrated intensity maps of HC^^N (4-3) (left panel) and SO (87 - Vy) (right panel) 
at the southern core. To avoid the influence of outflow motions, both the maps are integrated from 
2 km s"i to 8 km s"^ The contour levels are (a) -1.2 (60-), 1.2, 2.4, 4.2, 6.6, 9.6, 13.2, 17.4 
Jy beam"' -km s"' for HC'^N (4-3), (b) -1.2 (6(r), 1.2, 2.4, 4.2, 6.6, 9.6, 13.2, 17.4, 22.2, 27.6, 33.6 
Jy beam-i-km s-^ for SO (87 - 7?) 
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Fig. 10. — P-V diagrams and integrated intensity maps of SO (87 - 7?) (upper panels), and HC^^N 
(4-3) (lower panels) at the southern core. The P-V diagrams are cut along N-S direction. The 
vertical solid line in P-V diagrams labels the systemic velocity. The dashed and solid contours 
in the right panels show the red- and blue-shifted emission, respectively. The integral velocity 
intervals are marked by thick dashed lines in the P-V diagrams. For both SO (87 - Vy) and HC^^N 
(4-3), the blue-shifted emission is integrated from -4 km s"^ to km s ^ while the red-shifted 
emission from 10 km s'^ to 14 km s"^ in the integrated intensity maps, (a) P-V diagram of SO 
(87 - 77). The contours are from 0.6 (3cr) in steps of 0.6 Jy beam ^ (b) P-V diagram of HC^^N 
(4-3). The contours are from 0.6 (3cr) in steps of 0.4 Jy beam"^ (c) Integrated intensity maps of 
SO (87 - 77) at line wings. The contours are from 1 (Scr) in steps of 1 Jy beam"' -km s~' for both 
red- and blue-shifted emission, (d) Integrated intensity maps of HC'^N (4-3) at red wing. The 
contours are 0.6 (3cr), 1.2, 2, 3 Jy beam"' -km s"'. 
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Fig. 11. — CS (7-6) channel maps at the southern core, which is smoothed to a velocity resolution 
of 3 km s-\ The contours are -0.6 (3cr), 0.6, 1.2, 2.4, 4.8, 7.2, 9.6 Jy heanrK 
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Fig. 12. — HCN (4-3) channel maps at the southern core, which is smoothed to a velocity resolution 
of 4 km s-\ The contours are -0.6 (3cr), 0.6, 1.2, 2.4, 3.6, 4.8, 7.2 Jy beam'^ 
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Fig. 13. — P-V diagrams and integrated intensity maps of HCN (4-3) (upper panels), and CS (7-6) 
(lower panels) at the southern core. The P-V diagrams are cut along N-S direction. The vertical 
solid line in P-V diagrams labels the systemic velocity. The dashed and solid contours in the right 
panels show the red- and blue-shifted emission, respectively. The blue- and red-shifted emission 
in the integrated maps are integrated from -12 km s"^ to -5 km s ' and 15 km s"' to 22 km s \ 
respectively in (b) and (e) panels, (a) P-V diagram of HCN (4-3). The contours are from 0.9 (3cr) 
in steps of 1.2 Jy beam"^ (b) Integrated intensity maps of HCN (4-3) at line wings. The contours 
are 1.5 (5cr), 4.5, 7.5, 10.5 Jy beam~^-km s"^ (c) The integrated intensity maps of HCN (4-3) at 
extremely high velocities. The blue- and red-shifted emission in the integrated maps are integrated 
from -20 km s"^ to -13 km s"^ and 23 km s~^ to 39 km s \ respectively. The contours are from 1.5 
(5cr) in steps of 3 Jy beam"' -km s"' for both blue- and red-shifted emission, (d) P-V diagram of 
CS (7-6). The contours are from 0.9 (3cr) in steps of 1.2 Jy beam ^ (e) Integrated intensity maps 
of CS (7-6) at line wings. The contours are 1.5 (5cr), 4.5, 7.5, 10.5 Jy beam"' -km s"' 
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Fig. 14. — Population diagrams of H2CS towards four cm/mm cores. The names of the cores are 
labeled on the upper-right corner of each panel. Open circles in blue represent the observed data. 
The vertical bars present 3cr errors of ln(N„/g„) due to the uncertainties of integrated intensities. 
The solid line shows the linear least-squares fitting using the Rotational Temperature Diagram 
method. Crosses in red mark the weighted mean results from Population Diagram analysis. The 
inferred parameters from the Population Diagram analysis are presented on the upper-right comers 
of each panel. 
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Fig. 15. — Abundance ratios of [CS/SO] (left) and [CS/HCN] (right) versus flow velocity along the 
redshifted lobe. We range the excitation temperature from 10 K to 50 K to derive the abundance 
ratios of [CS/SO]. The excitation temperature in calculation of abundance ratios of [CS/HCN] is 
assumed to be 30 K. The solid line in the right panel is the linear least-squares fitting, and the 
fitting results are presented in the upper-right comer. 
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Fig. 16. — Mass- Velocity relationships for the outflow lobes. Left : blueshifted lobe; right: red- 
shifted lobe. The solid lines in both panels show the power law fit towards all the data. The 
dashed line in the right panel shows the power law fit towards the HCN and SO data up to 
V/tow = 25 km s~^ The fitting results are presented in the lower-left comers. 
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Table 1. Parameters of 860 fxm continuum emission 





R.A. 


Decl. 




Deconvolution sizes 


^peak 


S,, 


T/ 




Mass 




Name 


(J2000) 


(J2000) 




(" X ") 


(Jy beam" ' ) 


(Jy) 


(K) 




(Mo) 


(lO^^* cm-2) 


Northern core 


18:06:14.447 


-20:31:28.253 




Point source 


0.20±0.02 


0.26 


50 


1.5 


13 




Middle core 


18:06:14.668 


-20:31:31.830 


1.48" 


X 1.29" (P.A.=-37.8°) 


0.76±0.04 


1.07 


92 


1.2 


30 


1.2 


Southern core 


18:06:14.889 


-20:31:40.149 


4.71" 


X 1.26" (P.A.=-20.2°) 


0.95±0.12 


2.52 


51 


0.8 


165 


2.1 



"The dust temperature is assumed to be the same as the rotational temperature of H2CS transitions 
''The opacity index yS is obtained from Su et al. 2005 ) 



Table 2. Observed parameters of the lines 



Molecule 


Transition 


Frequency 


E„ 


rms 










Intensity'' 






fwhm'' 








(GHz) 


(K) 


(Jy beam"^) 




(km 








(Jy beam"' ) 






(km 
















D 


E 


F 


G 


D 


E 


F 


G 


D 


E 


F 


G 


H2CS 


10o,io-9o,9 


342.946 


90.6 


0.3 


5.5±0.2 


2.7±0.2 


5.9±0.2 


6.0±0.5 


1.9±0.3 


1.4±0.2 


1.7±0.1 


0.9±0.2 


2.7±0.4 


3.1±0.5 


4.3±0.6 


3.8±1.1 




102,9-92,8" 


343.322 


143.3 


0.3 


4.0±0.5 


2.2±1.1 




4.9±0.3 


0.9±0.2 


1.4±0.3 




1.1±0.3 


4.0±0.6 


4.4±1.6 




1.9±0.6 




102,8-92,7 


343.813 


143.3 


0.3 


6.0±0.3 


2.5±0.4 


4.9±0.5 




1.3±0.2 


0.9±0.1 


0.8±0.2 




3.0±0.7 


5.5±0.9 


3.9±1.3 






103,8-93,7" 


343.410 


209.1 


0.3 


6.2±0.8 


2.4±0.7 


5.0±0.3 




1.1±0.7 


l.liO.l 


2.1±0.3 




3.3±0.7 


4.1 ±0.5 


2.7±0.6 






103,7-93,6' 


343.414 


209.1 


0.3 


5.8±2.8 


2.6±0.2 


5.5±0.3 


5.8±0.2 


0.7±0.2 


2.0±0.2 


2.4±0.3 


1.5±0.3 


5.9±0.8 


4.0±0.5 


2.9±0.6 


2.3±0.5 




105.6/5-95,5/4^ 


343.203 


419.2 


0.2 




1.8±0.4 


5.9±0.3 


5.9±0.2 




0.6±0.2 


0.7±0.2 


0.9±0.2 




3.2±1.0 


2.4±0.8 


1.5±0.4 


SO 


83-77 


344.311 


87.5 


0.2 


4.9±0.1 


2.2±0.1 


5.0±0.1 


4.4±0.1 


2.9±0.1 


4.0±0.1 


5.5±0.1 


3.5±0.1 


3.9±0.2 


5.3±0.2 


9.0±0.2 


8.2±0.3 


HC'^N v=0 


4-3 


344.200 


41.3 


0.2 


6.1 ±0.4 


2.6±0.1 


5.2±0.1 


4.8±0.3 


0.6±0.1 


2.1±0.1 


2.8±0.1 


l.liO.l 


3.6±1.0 


4.8±0.3 


8.0±0.3 


7.5±0.8 


cs 


7-6 


342.883 


65.8 


0.3 


























HCN v=0 


4-3 


354.505 


42.5 


0.3 



























"Not all the detected lines are listed in this table. The others will be presented in another paper. 

''The V/jr, Intensity and FWHM of each transition are derived from single gaussian fit towards the beam-averaged spectra. 
•'Blended with H^^co (5i,4-4i,4) at 343.325713 GHz. 
■^Blended with H2CS (103,7-93,6) 
^Blended with H2CS (103,8-93,7) 

^The two transitions of H2CS (105,6-95,5) and (105,6-95,6) have the same frequency, line strength and permanent dipole moment. Therefore they has same contributions to the observed 
line profile. 



Table 3. The physical parameters of H2CS transitions obtained with Rotational Temperature Diagram (RTD) method and 

Population Diagram (PD) analysis 



Core RTD PD 





Trot (K) 


Nto, (1015 cm-2) 


Tro, (K) 


Nto, (lO'S cm-2) 


f 








T 






(10o,io-9o,9) 


(102,9-92,8) 


(102,8-92,7) 


(103,8-93,7) 


(103,7-93,6) 


(105,6/5-95,6/4) 


D 


43±9 


3.8±2.9 


42±34 


4.2±2.9 


0.46±0.24 


6.4±4.4 


0.7±0.5 


0.9±0.6 


0.4±0.5 


0.2±0.3 




E 


83±21 


2.5±1.6 


92±74 


3.6±3.0 


0.26±0.23 


4.1±4.3 


0.7±0.7 


0.6±0.6 


0.7±0.8 


0.7±0.8 


O.liO.l 


F 


83±7 


2.6±0.6 


51 ±23 


4.0±2.9 


0.34±0.23 


3.9±3.1 




1.1±0.9 


1.1±1.2 


1.1±1.2 


O.OiO.l 


G 


91±17 


1.3±0.7 


105±37 


3.7±3.1 


0.12±0.18 


2.5±2.7 


2.4±2.3 






2.4±2.1 


0.3±0.3 



^The rotational temperature and total column density of H2CS transitions derived from RTD analysis are presented in the second and third columns, while those derived from PD 
analysis are shown in the forth and fifth columns. The sixth column gives the filling factor of each source inferred from PD analysis. The last six columns exhibit the optical depth of each 
transition using PD analysis. 
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Table 4. Outflow parameters of the southern core 



Molecule Velocity interval M P E 

(km s-i) (Mo) (Mo • km s-^) (lO'^^erg) 



Component Blue Red Blue Red Blue red Blue Red 



SO [-4,0] [10,14] 13 13 86 82 5.8 5.4 

CS [-12,-5] [15,22] 3.7 5.5 47 68 6.0 8.7 

HCN [-20,-5] [15,39] 5.4 17.6 85 294 14.1 54.6 



